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A semi-classical Monte Carlo model for studying three-dimensional carrier dynamics in photocon- 
ductive switches is presented. The model was used to simulate the process of photo-excitation in 
GaAs-based photoconductive antennas illuminated with pulses typical of mode-locked Ti:Sapphire 
lasers. We analysed the power and frequency bandwidth of THz radiation emitted from these de- 
vices as a function of bias voltage, pump pulse duration and pump pulse location. We show that the 
mechanisms limiting the THz power emitted from photoconductive switches fall into two regimes: 
when illuminated with short duration (< 40 fs) laser pulses the energy distribution of the Gaus- 
sian pulses constrains the emitted power, while for long (> 40 fs) pulses, screening is the primary 
power-limiting mechanism. A discussion of the dynamics of bias field screening in the gap region 
is presented. The emitted terahertz power was found to be enhanced when the exciting laser pulse 
was in close proximity to the anode of the photoconductive emitter, in agreement with experimental 
results. We show that this enhancement arises from the electric field distribution within the emitter 
combined with a difference in the mobilities of electrons and holes. 

PACS numbers: 42.72.Ai, 73.20.Mf, 78.20.Bh, 78.47.+P 
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I. INTRODUCTION 

Over the last few years interest in using far-infrared 
radiation as a non-contact probe of the properties of 
various materials such as semiconductors [ll H, su- 
perconductors 1^ 1^ and biological tissue [3, 13 has 
grown considerably. This has been as a result of improve- 
ments in sources ]^ 0| and detectors of terahertz 
(THz) frequency radiation, and the development of novel 
spectroscopic methods such as time domain spectroscopy 
[O. Il3l| . By recording the electric field of a single cycle 
of terahertz radiation as function of time, the technique 
of time domain spectroscopy is capable of extremely sen- 
sitive measurements, especially in the far-infrared band. 
This technique represents a particularly powerful tool for 
studying the dynamics of weakly correlated systems in 
condensed matter physics, including excitons and Cooper 
pairs, which have binding energies observable in the THz 
region of the electromagnetic spectrum 0, |^ . 

Terahertz time domain spectroscopy relies on the abil- 
ity to produce pulses of THz radiation so short that they 
contain only one or even half a cycle. Single-cycle THz 
pulses can be generated by shining ultrashort laser pulses 
on bulk semiconductors such as InAs, GaAs and GaP 
[loL 13 ■ These pulses contain a continuous and broad 
distribution of frequencies, which can be exploited to 
perform spectroscopic studies over a large portion of the 
spectrum. 

In bulk semiconductors single-cycle radiation can be 
emitted via three mechanisms: the formation of an clcc- 
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trical dipole due to the difference of mobilities between 
electrons and holes (the photo-Dember effect) op- 
tical rectification [l^, and the acceleration of charges 
under the surface depletion field of a semiconductor ic|. 
Photoconductive switches (PCSs) are an attractive al- 
ternative to surface field semiconductors emitters for the 
production of single-cycle THz radiation ^3 • These de- 
vices have been shown to produce broad bandwidth THz 
signals with high electric field intensities |3- In a PCS 
two parallel metallic contacts are deposited on the sur- 
face of a semiconductor, an external biasing voltage is 
applied between these contacts and an ultra-short laser 
pulse is used to generate photo-carriers in the gap be- 
tween the contacts. The acceleration of these carriers by 
the external electric field gives rise to an electromagnetic 
emission according to Maxwell's equations. 

An understanding of the carrier dynamics in photo- 
conductive switches is of major importance in order to 
improve the bandwidth and power of these devices. Sev- 
eral theoretical and experimental studies have been per- 
formed to understand different aspects of the process of 
THz emission by PCS. These works include the anal ysis 
of charge transport by the use of rate equations 
and Drude-Lorentz models 0, |23 . In addition, the ef- 
fect of dislocations on the emitted signal has been studied 
by the inclusion of reduced mobilities in rate equations 
[2ll | as well as experimentally measuring time domain 
THz traces and fitting them to simple charge transport 
models _22J. Also, TDS experiments using double pulses 
to excite the PCS [IJ have been performed to obtain in- 
formation about the photo-generated carrier screening of 
the bias voltage. 

These models correctly predict most of the character- 
istics of the observed THz signals, but have the disadvan- 



tage that parameters such as mobihty and carrier trap- 
ping time have to be adjusted phenomenologically to re- 
produce the details of the experimental measurements. 
While such models have the advantage that under cer- 
tain circumstances they can be analytically solved, this 
possibility is limited by factors such as the geometry of 
the domain in which the equations can be solved, and 
in some cases the analysis has to be restricted to one or 
two dimensions [TsL [23. [2^ . A Monte Carlo simulation 
for PCSs was recently published that models the dy- 
namics of a low density of electrons generated by small 
photon fluences. However, it neglects the contribution of 
holes to the current density (see Section |n} as well as 
the carrier-carrier scattering mechanism which has been 
found to dominate the total carrier scattering rate . 

In this work a three-dimensional Monte Carlo model is 
presented (Section ^ that provides a powerful tool for 
the analysis of the influence of parameters on carrier dy- 
namics and the evolution of variables in time. The effect 
of the applied bias voltage (Section ^Oj and pump pulse 
width f Section II V|l on the emitted signal is analyzed. A 
discussion of electric field screening is also presented (Sec- 
tion Ej). 

II. MONTE CARLO MODEL 

In this paper we extend the Monte Carlo model of 
Ref. To] to simulate emission from PCSs. The model uses 
the current density vector J as source term of Maxwell's 
equations. The electric field is calculated by the use of 
the far-field approximation 

^ dJ 

Ethz oc — . 1) 
at 

In order to obtain the current density vector, the simu- 
lation volume is divided into a three-dimensional rectan- 
gular grid, and the displacement of each carrier (electron 
or hole) is calculated over small time intervals, which are 
shorter than the scattering times. In each time inter- 
val and spatial box, the electric field is assumed to be 
constant and uniform and then the displacement of each 
carrier is calculated by assuming it follows the solution 
of 



with initial conditions given by its previous position and 
velocity, where is the acceleration of the carrier with 
charge q and effective mass m* in an electric field E. 

Although the simulation assumes parabolic band dis- 
persion, inclusion of both T and L valleys creates a more 
realistic carrier ensemble. A set of quantum mechani- 
cal scattering rates are calculated. Random numbers are 
used to determine if each particle is scattered in each time 
step, the scattering angle and the energy loss. The scat- 
tering mechanisms included are: LO phonon absorption 
and emission within both T and L valleys, TO phonon 
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FIG. 1: (Colour online) Diagram of the boundary conditions 
used. On the plane z — the boundary conditions are set as 
for y < lO^m, Vb for y > 20/im and Neumann boundary 
condition is assumed in the region lO/xm < y < 20^m as well 
as on all other boundaries. 



emission and absorption between T and L valleys, acous- 
tic phonon, impurity and carrier-carrier scattering. De- 
tails of the implementation of these scattering mecha- 
nisms are described in Ref. [Tol| . 

In each time interval the charge density is calculated 
for each of the boxes in the spatial grid. This is then 
used to find the electric potential (j) by numerically solv- 
ing Poisson's equation in three dimensions using a relax- 
ation method with Chebyshev acceleration. The bound- 
ary conditions for solving Poisson's equation are set as 
follows (Fig. the potential is fixed to zero on one 
third of the top surface {z — 0), corresponding to one 
of the contacts, on the second third of the surface Neu- 
mann boundary conditions are assumed, and the poten- 
tial is fixed to Vb (the bias voltage) in the final third. 
Neumann boundary conditions are also assumed on the 
other boundaries. The electric-field vector E = — V(/), 
is determined using the central difference method. Ini- 
tially the spatial distribution of extrinsic carriers is cal- 
culated, and the temperature distribution of these cold 
carriers is set according to Maxwell-Boltzman statistics. 
The simulation is run until charge distribution equilib- 
rium is reached, prior to the photo-excitation process. 

Photo-excited carriers are assumed to be generated ac- 
cording to the spatial and temporal intensity of a Gaus- 
sian laser pulse absorbed in the sample. Random num- 
bers are used to calculate the initial positions of each 
photo-generated electron-hole pair. The distribution in 
the z-direction is weighted exponentially as defined by 
the absorption coefficient of the light in the semiconduc- 
tor. The initial energy of the electron-hole pairs is the 
difference between the semiconductor band gap and the 
photon energy, this energy is partitioned between elec- 
trons and holes based on momentum conservation. Pho- 
tons are assumed to have the energy distribution of a 
Gaussian pulse (thus spectral width in the simulation is 
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calculated from the laser pulse duration parameter). The 
directions of the pairs of momentum vectors are initially 
randomized with equal probability across 47r steradians. 
The net current density vector is then calculated as 



E 



qiri 



(3) 



its derivative is calculated and with it the radiated elec- 
tric field using the far field approximation (Equation^. 

The parameters used in the simulation are taken from 
Ref. uJ|. These parameters are typical of insulating 
GaAs (an n-type sample with very low donor density 
of ri = 2 X lO^^cm"'^ was chosen). A 4pJ laser pulse 
with 1.5 /im standard deviation beam- waist was simu- 
lated, where the centre wavelength of the laser pulse was 
chosen to be 800nm (1.55 eV). The domain of the sim- 
ulation was 30/xm x 30/xm x 5^m x 2ps, divided in to a 
space-grid of 64 x 64 x 32 and time intervals of 2 fs. The 
distance between the contacts of the simulated photocon- 
ductive switch was 10 /xm. One million pairs of pseudo- 
particles were used to represent extrinsic and intrinsic 
electrons and holes. 



III. 



BIAS VOLTAGE EFFECTS 



Using the model described in the previous section we 
performed a study of the effect of applied bias voltage 
on emitted THz power from GaAs PCSs. Since mode- 
locked TiiSapphire lasers are most commonly used in ex- 
periments we chose pulse parameters typical of this type 
of laser. The simulation results presented in Fig.[2Ji show 
the effect of altering the bias across a PCS. In this set 
of simulations the laser pulse duration was set to 10 fs, 
which is a pulse duration easily achieved using chirped 
mirror TiiSapphire oscillators. The THz pulse shortens 
at higher bias voltages as seen in Fig. O in which the 
FWHM of the pulses at different voltages are plotted. 
This is because at higher applied voltages a greater num- 
ber of electrons are transferred into the L- valley, resulting 
in a higher TO-phonon emission and absorption rate. As 
a consequence the average effective mass of the electron 
population increases (see inset in Fig. O, reducing the 
net electron mobility, thereby shortening the pulse. This 
effect has been observed in the experiments of Liu et al. 
[2^. In Fig. |2h the Fourier transforms of Fig. |2K are 
shown. These spectra broaden as the voltage increases, 
produced by the shortening of the THz transient. 

The power emitted increases with the bias voltage, as 
seen in Fig. |31 and in agreement with the experimen- 
tal measurements of This power enhancement is 
expected in a simple macroscopic picture. Taking the 
current vector J = crEfc, with a being the switch's con- 
ductivity and Eb the external bias field, assuming that 
Eft remain constant we get Ethz oc aEf,. If the power is 
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FIG. 2: (Colour online) Graphs of the THz waveform (upper 
panel) for lOfs pump pulses at different voltages, and their 
Fourier transforms (lower panel). The width of the pulse de- 
creases with voltage and this effect gives as a result an increase 
of the bandwidth. 



given by 
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and since |Eh| oc Vb we end up with P (x Vb ^, which is 
consistent with the quadratic-like behaviour of power in 
Fig. 121 Note that the assumption of E^ being constant, 
while a good first approximation, it is not strictly true 
as shown by the further analysis of the full microscopic 
simulation fSection lV)l . 



IV. PULSE WIDTH EFFECTS 

The excitation pulse duration is one of the main vari- 
ables affecting the temporal width (and hence the fre- 
quency bandwidth) of an emitted THz transient from a 
photoconductive switch. This is primarily because the 
pulse duration determines the time interval over which 
photo-carriers are generated, and therefore the current 
rise time in a PCS. Here we present a set of simulations 
where the excitation pulse duration is varied over a range 
achievable using TiiSapphire lasers (6 to 120 fs). 

In Fig.^the electric fields for various laser pulse dura- 
tions are shown for a bias voltage of 30 V. It can be seen 
that the THz transient becomes narrower as the pump 
pulse gets shorter (Fig.^^). In Fig.^ the spectra of the 
time-domain pulses are shown. While the distributions 
for long pump pulses (~ 120 fs) show a sudden drop in 
spectral components over 4 THz, distributions for short 
pulses (~ 10 fs) contain components over 8 THz. The 
distribution width at -lOdB of the spectra is used as a 
measure of the bandwidth of signal generated in the PCS 
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FIG. 3: (Colour online) Graph of full widths at half maxi- 
mum (circles) for 10 fs pump pulses as a function of the bias 
voltage showing the narrowing of the THz transient with volt- 
age. The power of the signal emitted (triangles) is also shown 
as function of applied voltage. The inset shows the evolution 
of the effective mass of electrons as a function of time for 0, 
10, 20 and 30V of bias voltage. 
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ingly a strong dependence can be seen, with the power 
dropping for both over and under 40 fs pulses. The drop 
for under 40 fs pulses can be explained macro scopically 
by calculating the current density vector, given by 



(5) 



Taking the number densities n and p to be proportional 
to the total number of carriers Nf, and iV/i, extracted from 
the simulation output, and calculating the mobilities fie 
and from the average effective mass, which is 
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for electrons and rrih for holes, the current vector can be 
estimated, and with it the emitted power, Nr and 
were again extracted from the simulation. In Fig. I^the 
dotted line shows the expected power using this method; 
the drop in power in the short pulse region approximately 
reproduces that of the simulation. It can therefore be 
attributed to the broadening of the energy distribution 
of the pump pulse as its duration gets shorter, having as 
a consequence two main effects. Firstly, the lower energy 
tail of the distribution falls below the band-gap energy 
(see inset in Fig. inj, resulting in a smaller density of 
photo-carriers generated. Secondly the high energy tail 
of the pump pulse distribution contributes to the power 
reduction, since a larger number of electrons are directly 
injected into the L valleys, resulting in a larger electronic 
effective mass and thus a lower mobility. 

In order to explain the drop in power for long pulses, 
it is necessary to take into account the actual electric 
field acting on the photo-generated carriers. This can be 
done by replacing the electric field in Equation|Slby an 
effective electric field obtained by averaging the simulated 
electric field over the size of the laser spot (^ ±0.8 /im 
from the center of the PCS's gap), and over ~ 100 fs (the 
screening time, see Section^. The dashed line in Fig.jS] 
shows the power calculated using this effective electric 
field. The good agreement between the simulations and 
this curve leads to the conclusion that the drop in power 
for long pulses can be attributed to the screening of the 
bias field. 



FIG. 4: (Colour online) Graphs of THz waveforms (upper 
panel) for different widths of the pump pulse durations. The 
THz pulse sharpens as the pump pulse width becomes smaller. 
The spectra of the previous signals is shown in the lower panel, 
there is a significant effect on the high frequency components. 

in Fig.[Sl This plot shows how the bandwidth of the sig- 
nal increases as the pump pulse shortens. 

A priori the power is not expected to depend signif- 
icantly on the exciting pulse width, as the number of 
photons per pulse and the applied bias voltage are iden- 
tical in the simulation. In Fig. [51 the power is plotted 
(triangles) as function of the pump pulse width. Surpris- 



V. SCREENING 

The power and bandwidth of radiation emitted by a 
PCS is influenced strongly by electric field screening ef- 
fects. Therefore it is not surprising that space-charge 
screening has been a major focus of models describing 
PCS emission, m 

Even prior to photoexcitation of a PCS, the steady- 
state distribution of extrinsic carriers around the elec- 
trodes must be considered in order to determine the com- 
plete electric vector-field within the active region of the 
PCS. In the present study this steady-state space-charge 
field is automatically calculated (before the optical exci- 
tation process is simulated) by iterating the simulation 
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FIG. 5: (Colour online) Graph of the spectral width at -lOdB 
as a function of pump pulse width (circles) at a bias voltage 
of 30 V. The width of the distribution increases as the pump 
pulse decreases caused by carriers being photo-generated in 
a shorter period of time. The power of the signal emitted is 
also plotted (triangles) as function of the pump pulse dura- 
tion, the dotted curve shows the power calculated taking into 
account the number of photo-generated carriers and their ef- 
fective mass, the dashed line includes additionally the effective 
bias electric field. The inset is a schematic with two curves 
(not drawn to scale) of the energy distributions of a short and 
a long pulse: A) are the photons with not enough energy to 
transfer carriers into the conduction band, B) is the band- 
gap energy, C) is the band-gap energy plus the L- valley offset 
and D) are the photons with enough energy to inject carriers 
directly into the L valley 
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FIG. 6; (Colour online) Graph of electric potential as func- 
tion of X and t. The surface show the screening of the electric 
field in the gap region after the arrival of the pump pulse (at 
t = 0) as the separation of electrons and holes forms a dipole. 



over enough time-steps until an equilibrium electric po- 
tential is reached. The form of the screened electric field 
between the electrodes for an n-type GaAs PCS may be 
inferred from the gradient of the t — — 1 ps trace in Fig.El 
It can be seen that the extrinsic electron distribution is 
such that the anode and cathode are screened to similar 
extents. 

Here we focus on the role of space-charge screening dur- 
ing and after the arrival of the pump pulse. On photo- 
excitation electron-hole pairs are generated within the 
high electric field gap region of a PCS. Charge separation 
produces a dipole, resulting in the emission of THz radi- 
ation. However, the resulting space-charge electric field 
opposes the applied field, thereby reducing the force ex- 
perienced by carriers. It is expected that this effect would 
be particularly significant for carriers photo-injected in 
the center of the PCS gap at later times. 

In Fig. Elthe simulated potential of the PCS is plotted 
(for z = 0) as a function of x-position and time. Before 
the arrival of the pulse (at ps) the gradient of the po- 
tential shows a continuous and positive x-component in 
the region of the gap between the electrodes. After the 
injection of photo-carriers the potential considerably flat- 
tens in the central part of the gap region, corresponding 
to screening of the electric field. 

It can be seen from Fig. that an abrupt change in 
the potential is present near the edges of the electrodes, 
which increases as a result of the photo-carrier-induced 
screening. Therefore, the gradient of the potential will 
have a large x-component, implying a high electric field 
zone. If the photo-carriers are generated near the elec- 
trodes an enhancement of the emitted signal may be ex- 
pected. 

Experimental studies have shown that the power emit- 
ted by a PCS may be enhanced by illuminating it asym- 
metrically. The enhancement is greatest if the region 
illuminated with the pump-pulse is close to the edge of 
the anode contact |3, 123, ISfl- This has important practi- 
cal implications for THz emitter design. 

Such "anode-enhanced" THz emission has been ex- 
plained previously b y a trap-enhancement of the electric 
field near the anode or by the presence of a high field 
region near the anode in the case of a Schot tky contact 
|31|. However, experimental measurements [30| suggest 
that the asymmetry between two Ohmic electrodes may 
be attributed to the difference of mobilities between elec- 
trons and holes. 

As discussed above, the simulation treats the case of 
PCS electrodes that are each screened to a similar ex- 
tent before photoexcitation. In reality, depending on the 
substrate and type of contacts, some asymmetry in this 
screening may exist. However here we show that the pres- 
ence of contact asymmetry is not a necessary condition 
in order to observe "anode-enhanced" THz emission. 

This effect is most easily illustrated by plotting the 
movement of charge over the duration of the simulation. 
The contour plots in Fig. IHa) and (b) show the change 
in charge density, Ap = p(300fs) — p(Ofs), as a function 
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of x-y position of the 10 /im gap PCS. Plots (a) and 
(b) are extracted from a simulation where the PCS was 
illuminated near the cathode and anode respectively. It 
can be seen that the displaced charge when the PCS is 
illuminated near the anode is considerable larger. 

The enhancement effect can be understood in terms of 
the higher effective mobility of the photo-injected elec- 
tron population over the hole population. When pho- 
toexcitation occurs close to the anode electrons (holes) 
migrate into the high (low) field region (see Fig. inj. 
Therefore, in this case it is the more mobile charges 
that are affected by the higher electric fields. Thus 
when electron-hole pairs are generated near the anode 
the higher electron mobility results in an increase in the 
effective current density, and as a consequence a larger 
amplitude THz electric field is emitted. 

The above explanation is somewhat simplified, as un- 
der typical electric fields F ^ L inter-valley scattering 
will lead to a drop in electron mobility. However, as 
the simulation considers intervalley scattering, all results 
presented here include this additional effect. 

A set of simulations were performed where the laser 
pulse arrival position was varied across the 10 /xm gap of 
the PCS. A summary of these data is shown in Fig.Cfc), 
where the peak electric field is plotted as a function of 
the distance of the laser pulse from the center of the PCS 
gap. The peak THz field is enhanced by approximately 
40% as the pulse approaches the edge of the anode. The 
field drops again nearer the anode as a larger fraction of 
the laser beam-waist is incident outside the active region 
of the PCS. These results are consistent with previous 
experimental observations. 

The widths of the Fourier transforms at -10 dB of the 
signal at different positions are also plotted in Fig. [7fc). 
It is interesting to note that the bandwidth increases as 
the excitation nears the anode. This can be attributed to 
the fact that in that region the electric field is greater, re- 
sulting in a shorter THz transient, as discussed in Section 

Finally, it is useful to characterise the timescale of 
space-charge screening during the photo-excitation of a 
PCS, so that the influence of pump-pulse duration on the 
THz emission process can be analysed. We fitted an error 
function to the x = ±2.5 fim. position slices of the graph 
shown in Fig.|H|and chose to define the screening time as 
2 standard deviations of that error function. The screen- 
ing time for this simulation was found to be ^ 100 fs for 
all simulated pump pulses between 10 and 100 fs. This 
suggests that space charge screening is dominated by the 
intrinsic dynamics of carriers in the semiconductor and 
it is not significantly affected by the carrier generation 
time. 



VI. CONCLUSIONS 

A semi-classical Monte Carlo model that simulates the 
carrier dynamics in a photoconductive switch was pre- 
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FIG. 7: (Colour online) Simulation results showing "anode 
enhanced" THz emission, (a) and (b) are contour plots of 
the change in charge density, Ap, for simulations where the 
laser pulse arrived near the cathode [x — — 3/im) and anode 
{x — 3/im) respectively. Ap is the change in charge density 
between t = to t = 300 fs projected onto the x-y plane. Ap 
was averaged over all slices in the z-direction. A summary 
of 11 simulations is shown in (c) where the peak THz electric 
field (triangles) and FWHM (circles) are plotted as a function 
of laser spot position. 



sented. The THz transient has been shown to narrow 
with the applied voltage, as well as to be strongly depen- 
dent on the pulse width of the exciting laser. The emit- 
ted THz power was also studied, and was found to have 
a maximum for pump pulses of 40 fs duration. A drop 
in power for shorter pulses was explained by a broad- 
ening in the energy distribution of the exciting pulse, 
which starts to contain a significant number of photons 
with sub-band-gap and over L- valley gap energies. This 
reduces the quantity of free photo-carriers and increases 
their average effective mass. The reduction in THz power 
for longer excitation pulses is explained by a decrease in 
the bias field as a result of screening, since the photo- 
carriers are injected over a longer period of time. 

The screening of the electric bias field was studied. It 
was found that the screening time is independent of the 
pump pulses (for pulses under 100 fs), implying that the 
screening time is mainly dependent on the carrier dynam- 
ics and not on their injection time on this timescale. The 
electric field was found to be screened over the gap be- 
tween the electrodes of the PCS, but not near the edges 
of the electrodes. This fact combined with the difference 
of mobilities of electrons and holes produces an increase 
in the peak THz electric field and a broadening in band- 
width when carriers are photo-excited near the anode 
edge. 

In the future this model promises to be useful for test- 
ing new PCS electrode designs and substrate materi- 
als (for example low temperature grown GaAs or ion- 
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implanted semiconductors) . For better correlation with wishes to thank CONACyT (Mexico) for a scholarship, 
experimental results it would also be informative to cou- 
ple a fidl three-dimensional electromagnetic simulation 
with the existing carrier dynamics simulation. 
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